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Motivated by recent experiments on cold atomic gases in ultra high finesse optical cavities, we 
consider the two-band Bose-Hubbard model coupled to quantum light. Photoexcitation promotes 
carriers between the bands and we study the interplay between Mott insulating behavior and super- 
fluidity. The model displays aU(l)xU(l) symmetry which supports the coexistence of Mott insulat- 
ing and superfluid phases, and yields a rich phase diagram with multicritical points. This symmetry 
is shared by several other problems of current experimental interest, including two-component Bose 
gases in optical lattices, and the bosonic BEC-BCS crossover for atom-molecule mixtures induced 
by a Feshbach resonance. We corroborate our findings by numerical simulations. 
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Introduction. — The spectacular advances in cold 
atomic gases have led to landmark experiments in 
strongly correlated systems. With the observation of the 
superfluid-Mott insulator transition in ^^Rb [ij, and the 
BEC-BCS crossover in 0], attention is now being 
directed towards multicomponent gases. Whether they 
be distinct atoms or internal states, such systems bring 
"isospin" degrees of freedom. They offer the fascinating 
prospect to realize novel phases, and to study quantum 
magnetism, Mott physics and superfluidity [3|. 

More recently, significant experimental progress has 
been made in combining cavity quantum electrodynam- 
ics (cavity QED) with ultracold gases. Strong atom-field 
coupling has been achieved using ultra high finesse op- 
tical cavities and with optical fibres on atom chips 
[^. These experiments open an exciting new chapter in 
coherent matter-light interaction, and have already led 
to pioneering studies of condensate dynamics Q. The 
light field serves not only as a probe of the many-body 
system, but may also support interesting cavity mediated 
phenomena and phases. This dual role has been exploited 
in studies ofpolariton condensates in semiconductor mi- 
eroeavities |7| . It is reinforced by ground breaking cavity 
QED experiments using superconducting qubits in mi- 
crowave resonators [1]. This has led to solid state mea- 
surements of the collective states of the Dicke model 0] , 
and remarkable observations of the Lamb shift ■ 

In this work we examine the impact of cavity radiation 
on the Bose-Hubbard model. We focus on a two-band 
model in which photons induce transitions between two 
internal states or Bloch bands. This is a natural general- 
ization of the much studied two-level systems coupled to 
radiation, and may serve as a useful paradigm in other 
contexts. The new ingredients are that the bosonic carri- 
ers may form a Mott insulator, or indeed condense. The 
primary question is whether a novel Mott insulating state 
can survive, which supports a condensate of photoexci- 
tations or mobile defects. In analogy with zero point 
motion in Helium [lH, this may be viewed as a form 



of supersolid in which fluctuations of the photon field 
induce defects. Whilst this question has its origin in po- 
lariton condensates in fermionic band insulators [l^ , the 
present problem is rather different. Since the integrity 
of the Mott state is tied to the interactions, a priori it 
is unclear that it survives the effects of itinerancy and 
photoexcitation. Nonetheless, the outcome is afhrmative, 
and the model displays both this novel phase and a rich 
phase diagram. Related phases were recently observed in 
simulations of other two-component models [l3. 14|. 

The Model. — Let us consider a two-band Bose- 
Hubbard model coupled to the light field of an optical 
cavity within the rotating wave approximation 



Jaajaj + Loip^ip + g ^ (bloiij + h.i 



(1) 



{ij)a 



where a 



a,b are two bands of bosons with [aijCcj] 



Sij . These might be states of different orbital or spin an- 
gular momentum. Here, Ca, effects the band splitting, Ua 
and V are interactions, Jq,, are nearest-neighbor hopping 
parameters, and uj is the frequency of the mode, tp. We 
consider just a single mode, which couples uniformly to 
the bands. The coupling, g, is the strength of the matter- 
light interaction. In view of the box normalization of the 
photon, we denote g = g/^/N, where N is the number of 
lattice sites. It is readily seen that Ni = ^'^d 



l)/2, commute with Hq. These 



are the total number of atomic carriers, and photoexci- 
tations (or polaritons) respectively. These conservation 
laws reflect the global U(l) x U(l) symmetry of Hq, such 
that a e'^a, h e'^^b, -0 e-'^^'f^, where i9, ip are 
arbitrary. This symmetry will have a direct manifesta- 
tion in the phase diagram, and suggests implications for 
other multicomponent problems. We begin by assuming 
that a are strongly interacting hardcore bosons, and that 
b are dilute so that we may neglect their interactions. 
Zero Hopping Limit. — To gain insig ht into (P) we ex- 



2 



amine the zero hopping limit. This will anchor the phase 
diagram to an exactly solvable many body limit. The 
photon couples all the sites, and in the thermodynamic 

limit is described by a coherent state, I7) = e~~ e^^ |0), 
with mean occupation (ip'^ip) = 7^. We may replace the 
grand canonical Hamiltonian, H = Hq — fJ.iNi — /i2-^2, 
by an effective single site problem, (7|-ff I7) = J2i '^i- 



+ gj{b^a + a''b), 



(2) 



and we drop the offset, — /i2/2. We define ia = — 
A^i + A*2/2, Cb = e& - /ii - /i2/2, (Ij = cj - /i2, and the 
mean photon occupation per site, 7^ = 7^/iV. The effec- 
tive Hamiltonian ([2]) describes a single two-level system 
coupled to an effective "radiation field" of &-particles, or 
the Jaynes-Cummings model; for N two-level systems 
this is known as the Dicke or Tavis-Cummings model, 
and is integrable [H, These paradigmatic mod- 

els are well known in both atomic physics and quan- 
tum optics, and describe localized excitons coupled to 
light [131 . The eigenstates of ^ are superpositions in 
the upper and lower bands (that we denote as |0,n) 
and \l,n — 1)) with total occupancy n. The lowest 
energy is £',7 = ui-=f^ + net - <io/2 - ^12)^/4 + g^j'^n, 
where Cbo = et ~ ^a- Minimizing on 7 gives a self- 
consistency equation for the photon, and the result- 
ing eigenstates yield the zero hopping phase diagram in 
Fig. [TJ In the thermodynamic limit described here, only 
the lowest Mott state, with Ua + ni, = 1, survives; for 
Ml 5^ £b ~ /^2/2 — g'^/4:Uj it is favorable to macroscopi- 
cally populate the upper band. Within this stable Mott 
phase the variation yields either 7 = 0, corresponding to 
zero photon occupancy, or 7^ = {g^ — ^^)/4a)^g^, where 
(jc = -n/Swq; the former occurs for g < gc and the latter 
for 9 ^ 9c- In fact, this onset corresponds to the super- 
radiance transition in the Dicke model 3, l3 
since Ha + ni, = I = 2S in the lowest lobe, 



. Indeed, 
one may 

construct the Dicke model directly from ([1]) by using 
a spin 5 = 1/2 Schwinger boson representation, where 

S+ = b^a, S- = a'^b, = {m - na)/2: 



(3) 



This describes N two-level systems (or spins) coupled 
to photons, and may be treated using collective spins, 
J = Si. This yields a large effective spin, which 
may be treated semiclassically as — > 00. The on- 
set of the photon is accompanied by a magnetization, 
Ai = {J^)/N, which also serves as an order parameter 
for this continuous transition: A4 = —1/2, for g < gc, 
and M = —{gc/g)'^/2, for g > gc- This growth reflects 
the population imbalance, (rib) — (na), due to photoex- 
citations. The agreement between the variational and 
Dicke model results is a useful platform for departures. 

Variational Phase Diagram. — Having confirmed a zero 
hopping Mott phase, with Ua + Ub = 1, let us consider 
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FIG. 1: Zero hopping phase diagram in the large- C/a limit, 
with £(, = —ea = uj = g — 1, corresponding to tJ < wo- The 
vertical line, 3 = gc, is the superradiance transition in the 
Dicke model, and separates a Mott insulator with Ua + rib — 1 
and = (dark blue), from a superradiant Mott insula- 

tor with {tp^ij}) 7^ (light blue). Outside are the vacuum, and 
an unstable region corresponding to macroscopic population 
of the b states. Whilst the total density is fixed within both 
Mott phases, the individual a and b populations vary in the 
superradiant phase and may be viewed as isospin order. For 
uj > ujQ, the boundaries may cross and terminate the lobe. 



itinerancy and carrier superfluidity. Within this lowest 
lobe, we may consider hardcore a and b bosons [2^. A 
convenient approach is to augment the variational anal- 
ysis of Ref. y with a coherent state for the light field: 

|V) = I7) «) n [cos 61, (cos x»al + sinxi^l) 



4- sin 6*^ (cos ryj -I- sin r/i&|a|) |0) 



(4) 



where I7) is the coherent state introduced previously, and 
X; '7j 7 are to be determined. The first term in brackets 
describes the Mott state, and the second superfluidity. 
For 9 = this coincides with the variational approach 
for localized excitons coupled to light [l^l and reproduces 
the previous results for = 0. More generally, (g]) takes 
into account the effects of real hopping, involving site va- 
cancies and interspecies double occupation. It provides 
a starting point to identify the boundaries between the 
Mott and superfluid regions. We consider spatially uni- 
form phases, with energy density, £ = (V\H\V) /N: 

£ = {e+ - e_ cos 2x) cos^ 6 + (2e+ + V) sin^ 77 sin^ 9 

- ^ [Ja cos^ (x - v) + Jb sin^ (x + v)] sin^ 29 

+ 0)7^ -I- 57 cos^ 9 sin 2x, 

_ _ _ (5) 

where z is the coordination number and e± = (eb±eo)/2. 
Minimizing on 7 yields 7 = — g cos^ sin 2x/2a;, and 
one may eliminate this from £. Exploiting symmetries 
one may minimize £ over [0, 7r/2]. The order parame- 
ters, (a) = i sin20cos(x — ?y), {b) = 5 sin26'sin(x -I- ??), 
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FIG. 2: Variational phase diagram with Ja = Jb = J and 
ea — —1, eb = 1^ ~ V = g = 1, = 0. The phases 
are (i) a Mott insulator (dark blue), (ii) a superradiant Mott 
state supporting a condensate of photoexcitations (light blue), 
(iii) a superradiant superfluid (light red), and (iv) an a-type 
superfluid (dark red). The circles denote the transition to 
superfluidity as determined by 6, and the squares denote the 
onset of photons as determined by x- For these parameters, 
the transition from the Mott insulator to the superradiant 
Mott state occurs for = (3 - V5)/2 ^ 0.382. 
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and {ip^^p)/N = 7^, yield the phase diagram in Fig. [21 
where Ja = Jb = J- For the chosen parameters, we 
have four distinct phases; (i) a Mott state with (a) ~ 
(&) = (^/'^^/') = 0, (ii) a superradiant Mott state with 
(a) — (b) = and (ip'^ip) 7^ 0, (iii) a single component 
superfluid with (a) 7^ and (6) = (ip'^ip) — 0; (i^) 
a superradiant superfluid (a) 0, (b) ^ 0, (V'^V') 7^ 0- 
Indeed, the Hamiltonian displays a U(l) x U(l) symme- 
try and these may be broken independently. The phase 
diagram reflects this pattern of symmetry breaking. In 
particular, the superradiant Mott state corresponds to an 
unbroken U(l) in the matter sector (corresponding to a 
pinned density and phase fluctuations) but a broken U(l) 
(or phase coherent condensate) for photoexcitations. The 
expectation value of the bilinear^ (^^a) 7^ 0, corresponds 
to the onset of coherence in the Dicke model. This novel 
phase may be regarded as a form of supersolid. 

In the absence of competition from other phases, the 
transition between the non-superradiant insulator (Q = 
X = 7 = 0) and the a-type superfluid {0 ^ 0, x ^ V = 
7 = 0) occurs when Ca + zJ = 0. In Fig. [21 this is the line 
fi2 = 2(1 — z J). This crosses the superradiance boundary 
at a tetracritical point {zJ^^is'^) = ('^/2,2 — r), where 
r = (1 + -^5)72 is the Golden ratio. This follows from a 
Landau expansion of (O; eliminating 7, all the quadratic 
"mass" terms vanish. More generally, the phase diagram 
evolves with the parameters, and the a-type superfluid 
may be replaced by the proximate phases [l7[. 

Numerical Simulations. — To corroborate our findings, 
we perform exact diagonalization on a ID system of hard- 
core a and b bosons, with iV = 8 sites and periodic 
boundary conditions. The Hilbcrt space is truncated to 
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FIG. 3: Exact diagonalization for a ID system with 8 sites, 
= 16 photons, and the parameters of Fig. [2] We show 

(a) the total atom density and the Mott-superfluid transition, 

(b) the density of photons (reduced by a factor of two) and 
the onset of supperradiance, (c) the density of a-atoms, and 
(d) the density of fo-atoms. The dashed lines are a guide to 
the eye, and indicate the Mott-superfluid and superradiance 
transitions, as determined by hand from (a) and (b). Their 
intersection yields the location of the tetracritical point. 



a maximum number of photons, ~ 2N = 16. Fig. [3] 
shows the total atom, photon, a-atom and fe-atom den- 
sity. The dashed lines indicate the approximate locations 
of the Mott-superfluid (vertical line) and superradiance 
(horizontal line) transitions, as determined from panels 
(a) and (b). Although an accurate phase diagram for the 
thermodynamic limit is beyond the scope of this work, 
the features are in excellent agreement with Fig. [21 This 
parallels the success of mean field theory in other low- 
dimensional bosonic systems, and is remarkable given the 
enhanced role of fiuctuations. This may be assisted by 
the long range nature of the cavity photons. The super- 
radiance transition encompases the superfluid and Mott 
phases, and yields a tetracritical point; see (a) and (b). 
In addition, the region of a-density over-extends that of 
fe-dcnsity resulting in a pure a-type superfluid; see (c) 
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and (d). Our simulations suggest that this phase is sta- 
ble with increasing system size [17 1 . 

Discussion. — A feature not addressed by the present 
mean field theory, but captured in Fig. [31 is the disper- 
sion of the superradiance transition with J; in the Mott 
phase, 9 = 0, and J drops out of the variational energy 
One way to understand this is to recast the matter 
contribution as |Vm) = Oilcosxi + sin Xiblai)\fl), where 
I $7) = Yii al|0). This only accomodates Zoca/ particle-hole 
pairs. By analogy with the BCS approach to exciton in- 
sulators the Mott state may be refined and the J 
dependence restored by incorporating momentum space 
pairing 17| . This connection to the BEC-BCS crossover 
for bosons is reinforced by the Feshbach resonance 

0, 



21, 2 



problem studied in the absence of a lattice 
Performing a particle-hole transformation, the matter- 
light coupling reads ip^aibi. Aside from the global nature 
of the photon, this converts a and b into a "molecule" "0. 
At the outset there are eight phases corresponding to con- 
densation of (a), (b), (t/)). Of these, only five may survive; 
condensation of two variables provides an effective field 
(as dictated by the coupling) which induces condensation 
of the other. The band asymmetry, ea < Ch, reduces this 
to four, or less, depending on the parameters. In contrast 
to the single species mean field theory, this case supports 
an atomic superfluid, since condensation of one carrier no 
longer induces a field. Moreover, condensation may leave 
a U(l) symmetry intact, which allows the coexistence of 
Mott insulating and phase coherent behavior. 

In deriving Q and the phase diagram, we are primar- 
ily concerned with the matter-light coupling. As such, 
we incorporate V as in Ref. |3. This gives rise to the non- 
trivial phases in Fig. [2l However, as noted by Soyler et 
al [l3| . analogous phases may be stabilized in the two- 
component Bose-Hubbard model, without matter-light 
coupling, through a more sophisticated treatment of V it- 
self. Indeed, onsite repulsive interactions, Vnafii,, favor a 
particle of one species and a hole of the other on the same 
site. Treating this pairing in a BCS approach, one may 
replace n°n^ by |Aip + (Ai6jai+h.c.), where A.^ = (a|6i), 
is to be determined self-consistently. This field acts as 
a local "photon" , and a similar mean field phenomenol- 
ogy may ensue. Such pairing also occurs in fermionic 
models Although our discussion has focused on a 

single global photon, the symmetry analysis is more gen- 
eral. This is supported by studies of the two-band Bose- 
Hubbard model for equal fillings and commensurate den- 
sities |2J| . We shall provide details of the similarities and 



differences of this local problem in Ref. [17|. The classical 
limit may also be realized in optical superlattices, where 
giOibl is tunnelling between wells, a and b. 

Conclusions. — We have considered the impact of pho- 
toexcitations on the Bose-Hubbard model. The phase 
diagram supports a novel phase where photoexcitations 
condense on the background of a Mott insulator. We 



have performed numerical simulations, and highlighted 
connections to other problems of current interest. Di- 
rections for research include the impact of fluctuations 
and the nature of collective excitations. It would also 
be interesting to incorporate a finite photon wavevector. 
This may stabilize inhomogeneous phases and probe in- 
commensurate magnetism. Recent studies of Bose- Fermi 
mixtures also display a similar phenomenology, in which 
superfluidity is replaced by fermionic metalicity [l^ . 
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